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Abstract 

The rate of gravitational- wave energy loss from inspiralling binary systems of 

compact objects of arbitrary mass is derived through second post-Newtonian 

(2PN) order 0[(Gm/rc 2 ) 2 ] beyond the quadrupole approximation. The result 

has been derived by two independent calculations of the (source) multipole 

moments. The 2PN terms, and in particular the finite mass contribution 

therein (which cannot be obtained in perturbation calculations of black hole 

spacetimes), are shown to make a significant contribution to the accumulated 

phase of theoretical templates to be used in matched filtering of the data from 

future gravitational-wave detectors. 
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One of the most promising astrophysical sources of gravitational radiation for detection 
by large-scale laser-interferometer systems as the US LIGO or the French-Italian VIRGO 
projects |1| is the inspiralling compact binary. This is a binary system of neutron stars or 
black holes whose orbit is decaying toward a final coalescence under the dissipative effect of 
gravitational radiation reaction. For much of the evolution of such systems, the gravitational 
waveform signal is an accurately calculable "chirp" signal that sweeps in frequency through 
the detectors' sensitive bandwidth, typically between 10 Hz and 1000 Hz. Estimates of 
the rate of such inspiral events range from 3 to 100 per year, for signals detectable out to 
hundreds of Mpc by the advanced version of LIGO @. 

In addition to outright detection of the waves, it will be possible to determine parameters 
of the inspiralling systems, such as the masses and spins of the bodies to measure 

cosmological distances ||, to probe the non- linear regime of general relativity j7|, and to 
test alternative gravitational theories || . This is made possible by the technique of matched 
filtering of theoretical waveform templates, which depend on the source parameters, against 
the broad-band detector output ||. 

Roughly speaking, any effect that causes the template to differ from the actual signal by 
one cycle over the 500 to 16,000 accumulated cycles in the sensitive bandwidth will result in 
a substantial reduction in the signal-to-noise ratio. This necessitates knowing the prediction 
of general relativity for gravitational radiation damping, and its effect on the orbital phase, 
to substantially higher accuracy than that provided by the lowest-order quadrupole, or 
Newtonian approximation. If post-Newtonian corrections to the quadrupole formula scale 
as powers of v 2 ~ m/r (G = c = 1), then, say, for a double neutron-star inspiral in the 
LIGO/VIRGO bandwidth, with m/r typically around 10 -2 , corrections at least of order 
(m/r) 2 ~ 10 -4 will be needed in order to be accurate to one cycle out of the 16,000 cycles 
accumulated for this process. This corresponds to corrections of second post-Newtonian 
(2PN) order. 

Although numerous corrections to the quadrupole energy-loss formula have previously 
been calculated (for a summary, see fL0|), the 2PN contributions to the energy loss rate 
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for arbitrary masses have not been derived. This paper presents these contributions for the 
first time, discusses their significance for gravitational- wave data analysis, and outlines the 
derivation. 

The central result is this: through 2PN order, the rate of energy loss, dE/dt, from a 
binary system of two compact bodies of mass m% and m 2 , orbital separation r, and spins Si 
and S2, in a nearly circular orbit (apart from the adiabatic inspiral) is given by 
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where L is a unit vector directed along the orbital angular momentum, m = mi + rri2, 
i] = m 1 m 2 /m 2 , Xi = ^>i/ m \i X2 = ^2/ m 2 an< i denotes the sum over % = 1,2. The terms 
in square brackets in Eq. are respectively: at lowest order, Newtonian (quadrupole); at 
order m/r, 1PN [[ll]]; at order (m/r) 3//2 , the non-linear effect of "tails" (Air term) |]P2| , |l3"|| , 
and spin-orbit effects [ TA ,P^3 1 ; and at order (m/r) 2 , the 2PN terms (new with this paper), 



and spin-spin effects |i4|- |lq1 . 

For the special case of a test mass orbiting a massive black hole, perturbation theory 
has been used to derive an analogous analytic formula (apart from spin-spin effects) Ul2| , |17f . 
and for non-rotating holes, to extend the expansion through the equivalent of 4PN order 
|18| . The test-body 77 = limit of Eq. ([!]) agrees completely with these results to the 
corresponding order. 

The equations of motion for circular orbits, correct to 2PN order including spin effects, 
yield for the orbital angular velocity uo = v/r and the orbital energy E pT9| , [T^JT0[ 
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Combining Eqs. (0) and (0), one can express the rate of change u of the angular velocity as 
a function of lu, and get 
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where the spin-orbit {(3) and spin-spin (cr) parameters are given 

75t7)L • %j , and cr = (77/48) (—247%! ■ % 2 + 721L • Xi^ 1 ' ^2)- From that one calculates 
the accumulated number of gravitational- wave cycles M = J(f/f)df, where / = uj/tt is 
the frequency of the quadrupolar waves, in terms of the frequencies at which the signal 
enters and leaves the detectors' bandwidth. In order to avoid complications caused by spin- 
induced precessions of the orbital plane pO"|Jl5"]| , we assume that the spins are aligned parallel 
to the orbital angular momentum (in particular j3 and a remain constant). Using 10 Hz as 
the entering frequency of LIGO/VIRGO-type detectors, set by seismic noise, and, as the 
exit frequency, the smaller of either 1000 Hz (set by photon-shot noise) or the frequency 
corresponding to the innermost stable circular orbit and the onset of plunge [for small mass 
ratio, / max = l/(6 3 / 2 7rm) |[21|| 1, we find contributions to the total number of observed wave 
cycles from the various post-Newtonian terms listed in Table 1. 

Because % < 1 for black holes, and ^ 0.63 — 0.74 for neutron stars (depending on the 
equation of state, see ||22|), (3 and a are always less than ~ 9.4 and ~ 2.5, respectively. 
However, if we consider models for the past and future evolution of observed binary pulsar 
systems such as PSR 1534+12 and PSR 1913+16, we find (using a conservative upper limit 



for moments of inertia) that Xi 534+12 < 5.2 x 1(T 3 , Xi 913+16 < 6.5 x 1(T 3 , and we expect 
X2 ~ 7 x 10~ 4 . If such values are typical, both the spin-orbit and (a fortiori) the spin-spin 
terms will make negligible contributions to the accumulated phase. 

Table 1 demonstrates that the 2PN terms, and notably the finite-mass (//-dependent) 
contributions therein (which cannot be obtained by test-body approaches), make a significant 
contribution to the accumulated phase, and thus must be included in theoretical templates 
to be used in matched filtering. The additional question of how the presence of these terms 
will affect the accuracy of estimation of parameters in the templates can only be answered 
reliably using a full matched filter analysis @||. This is currently in progress ||23|| . 

The remainder of this paper outlines the derivations leading to this result. Two entirely 
independent calculations were carried out, using different approaches, one by BDI, using 



their previously developed generation formalism [24,25], the other by WW, using a formal 



slow-motion expansion originated by Epstein and Wagoner p6) . Details of these calculations 



will be published elsewhere p7j . 



Both approaches begin with Einstein's equations written in harmonic coordinates (see 
for definitions and notation). We define the field h a ^, measuring the deviation of the 
"gothic" metric from the Minkowski metric r] al3 = diag(— 1, 1, 1, 1): h al3 = \/—gg a ^ — r) a ^. 
(Greek indices range from to 3, while Latin range from 1 to 3). Imposing the harmonic 
coordinate condition dph a ^ = then leads to the field equations 

Uh aP = 167i(-g)T aP + A a(3 (h) = 16nr al3 , (4) 

where □ denotes the flat spacetime d'Alembertian operator, T af3 is the matter stress-energy 
tensor, and A a @ is an effective gravitational source containing the non-linearities of Einstein's 
equations. It is a series in powers of h a/3 and its derivatives; both quadratic and cubic 
nonlinearities in A a/3 play an essential role in our calculations. 

Post- Minkowski matching approach (BDI). This approach proceeds through several steps. 
The first consists of constructing an iterative solution of Eq. (|j) in an inner domain (or near 
zone) that includes the material source but whose radius is much less than a gravitational 



wavelength. Denning source densities a = T 00 + T kk , Oi = T 0t , <Jij = T % \ and the retarded 
potentials V = -47rD^V, Vi = -47rD^V i , and Wij = -47rD^ 1 [cr ij + (4Ti)~ 1 (d i Vd j V - 
hSijdkVdkV)), where n^ 1 denotes the usual flat spacetime retarded integral, one obtains the 
inner metric to some intermediate accuracy 0(6, 5, 6): h®® = —4V + 4(Wu — 2V 2 ) +0(6), 
= —4V; + 0(5), hfa = —4:Wij + 0(6), where 0(n) means a term of order e n in the post- 
Newtonian parameter e ~ v/c. From this, one constructs the inner metric with the higher 
accuracy 0(8, 7, 8) needed for subsequent matching as = □^ 1 [16vrr Q/3 (V r , W)}+0(8, 7, 8), 
where T a ^{V, W) denotes the right-hand-side of Eq. @j when retaining all the quadratic and 
cubic nonlinearities to the required post-Newtonian order in the near zone, and given as 
explicit combinations of derivatives of V, Vi and W^. The second step consists of constructing 
a generic solution of the vacuum Einstein equations (Eq. (£|) with T af3 = 0), in the form of a 
multipolar-post-Minkowskian expansion that is valid in an external domain which overlaps 
with the near zone and extends into the far wave zone. The construction of h a ^ in the 
external domain is done algorithmically as a functional of a set of parameters, called the 
"canonical" multipole moments M i± ...»,(£), ^...^(i) which are symmetric and trace-free (STF) 
Cartesian tensors. Schematically, h^ t = T a ^\M L , Sl] where L = i x ...%i and where the 
functional dependence includes a non-local time dependence on the past "history" of Mi(t) 
and Si{t). The third, "matching" step consists of requiring that the inner and external 
metrics be equivalent (modulo a coordinate transformation) in the overlap between the inner 
and the external domains. This requirement determines the relation between the canonical 
moments and the inner metric (itself expressed in terms of the source variables). Performing 
the matching through 2PN order p5[ thus determines M^(t) = /L[r a/3 ] + 0(5), Sl(£) = 
j L pfa/3] +0(4), where the "source" moments II and Jl are given by some mathematically 
well-defined (analytically continued) integrals of the quantity r a/3 (V, W) which appeared as 
source of h^. When computing the source moments we neglect all finite size effects, such 
as spin (which to 2PN accuracy can be added separately) and internal quadrupole effects. 
The final result for the 2PN quadrupole moment reads, in the case of a circular binary, 
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Iij = STFy 77m [ArV + BrVv 5 '] , 



(5) 



with A = 1 - i(m/r)(l + 39r/) - T ^(m/r) 2 (461 + 18395r/ + 24V), B = g(l - 3r/) + 
gig(m/r)(1607 — I68I77 + 229?7 2 ). The final step consists of computing from the external 
metric h^t the gravitational radiation emitted at infinity. This entails introducing a (non 
harmonic) "radiative" coordinate system X M = (T, X 1 ) adapted to the fall-off of the metric 
at future null infinity. The transverse-traceless (TT) asymptotic waveform hf? can be 
uniquely decomposed into two sets of STF "radiative" multipole moments Ul, Vl which are 
then computed as some non linear functionals of the canonical moments, and therefore of 
the source multipole moments. For instance, up to 0(5), 



which contains a non-local "tail" integral (in which b\ = be~ u ^ 12 where b is a freely specifiable 
parameter entering the coordinate transformation x^ — > : T — R = t— |x| — 2m ln(|x|/6)). 
The superscript (n) denotes n time derivatives. The energy loss is given by integrating the 
square of dhf T /dT over the sphere at infinity. At 2PN order this leads to 



Inserting the 2PN expression (|) of the source quadrupole into the radiative quadrupole (0), 
and using the previously derived 1PN expressions for the other multipole moments, we end 
up with the energy loss ([!]). 

Epstein-Wagoner approach (WW). This approach starts by considering 
_ n^ 1 (lQ7iT alS ) as a formal solution of Eq. (^) everywhere. One then expands the re- 
tarded integral d^ 1 to leading order in 1/R in the far zone, while the retardation is expanded 
in a slow-motion approximation. Using identities, such as = ^d 2 (r 00 x l x^) / 'dt 2 + spatial 
divergences, which result from dpT a " = (a consequence of the harmonic gauge condition), 
we express the radiative field as a sequence of Epstein- Wagoner multipole moments, 
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h %T = ~ E • • ■ ^^•"* m (< - R)tt , (8) 

at m =0 

where Jew ™ are integrals over space of moments of the source r Q/3 (e.g. J^w = 
/ T 00 x l x^d 3 x; see [26,28] for formulae). To sufficient accuracy for the radiative field, the 
2-index moment must be calculated to 2PN order, the 3 and 4-index moments to PN order, 
and the 5 and 6-index moments to Newtonian order. The moments of the compact-support 
source distribution (—g)T a/3 are straightforward. Contrary to what happens in the BDI 
calculation where the matching leads to mathematically well-defined formulas for the source 
multipole moments, the E W moments of the non-compact A a/3 source are given by formally 
divergent integrals. To deal with this difficulty we define a sphere of radius 1Z ^> A ~ r/e 
centered on the center of mass of the system, and integrate the non-compact moments 
within the sphere. Many integrations by parts are carried out to simplify the calculations, 
and the resulting surface terms are evaluated at 1Z and kept. The divergent terms are 
proportional to 1Z, and signal the failure of the slow-motion expansion procedure extended 
into the far zone. We discard the divergent terms. In order to compare directly with BDI, 
we transform the EW moments into STF moments using the projection integrals given by 
Thorne [p9[| . For the quadrupole moment, for instance, that transformation is given by 

jij _ Qr-r-p r r«? i 1 f-t 1 jijkk -, jkijk , jkkijs . 1 /r>o jijkkll qo jkijkll . n jkkijll <-> j-klklijxi 
1 — Oir ij^EW" 1 " 21v iiJ EW izJ EW "T^EW J "T 63^ ° EW OZJ EW + iui EW +zi EW )\- 

We find that those STF moments agree exactly with the "source" moments of BDI, e.g. 
Eq. (|5|). Note that the formal EW approach misses the tail effects in the waveform (see (|6|)). 
They must be added separately. 
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TABLES 

TABLE I. Contributions to the accumulated number N of gravitational-wave cycles in a 
LIGO/VIRGO-type detector. Frequency entering the bandwidth is 10 Hz (seismic limit); fre- 
quency leaving the detector is 1000 Hz for 2 neutron stars (photon shot noise), and ~ 360 Hz and 
~ 190 Hz for the two cases involving black-holes (innermost stable orbit). Spin parameters (3 and 
a are defined in the text. Numbers in parentheses indicate contribution of finite-mass (??) effects. 



2 x 1.4M© 1OM + 1.4M 2 x 10M ( 

Newtonian 16,050 3580 600 

First PN 439(104) 212(26) 59(14) 

Tail -208 -180 -51 

Spin-orbit 17(3 U(3 A/3 

Second PN 9(3) 10(2) 4(1) 

Spin-spin —la — 3<7 —a 



12 



